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Abstract. The energy carried in the gravitational wave signal from an eccentric bi- 
nary is spread across several harmonics of the orbital frequency. The inclusion of the 
harmonics in the analysis of the gravitational wave signal increases the signal-to-noise 
ratio of the detected signal for binaries whose fundamental frequency is below the 
galactic confusion-limited noise cut-off. This can allow for an improved angular reso- 
lution for sources whose orbital period is greater than 2000 s. Globular cluster sources 
includ possible binary black holes and neutron stars which may have high eccentricities. 
Cluster dynamics may also enhance the eccentricities of double white dwarf binaries 
and white dwarf-neutron star binaries over the galactic sources. Preliminary results of 
the expected signal-to-noise ratio for selected globular cluster binaries are presented. 



The space-based gravitational radiation detector LISA will be sensitive to fre- 
quencies in the millihertz range. A number of interesting astrophyisical phenomena 
will be likely sources in this frequency band, many of interest in cosmology and as 
fundamental tests of relativity. LISA will also be of use as an astronomical tool for 
observing compact binary systems in the galaxy and globular cluster system. In 
the globular cluster system, LISA can test dynamical evolution nodels by detecting 
binary populations which are predicted by simulations but which are exceedingly 
difficult to observe in electromagnetic radiation. One such class of objects is a pop- 
ulation of binary black holes with separations between 3 and 30 Rq and a thermal 
distribution of eccentrities [1,2]. A Monte Carlo analysis of the detection of such 
signals indicate that eccentricities above 0.8 can extend the detectable systems to 
orbital periods on the order of 1 day. 

The gravitational wave signal from an eccentric binary is given by the two po- 
larizations [3]: 
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in a coordinate system centered on the binary with the x-axis along the semi-major 
axis and the z-axis along the angular momentum vector of the system, with the 
orbital frequency of the binary given by /. The angles i9 and (/? are the usual polar 
coordinates giving the direction to the detector. The metric components, h^^^ and 

h^x±yj share a common amplitude which depends upon the masses and the orbital 
frequency. In general, they are combinations of Bessel functions. 

The signal detected by LISA is modified by the motion of LISA in its orbit about 
the sun. To describe the detected signal, we define several unit vectors. Let the 
orientation of the arms of the constellation of LISA spacecraft be given by the three 
vectors, ^i, £2, £3- The direction from the sun to the source is n. The orientation of 
the binary orbit is L (which points along the angular momentum) , d (which points 
along the semi- major axis), and b — L x d. The gravitational wave propagates 
along —n, and we define the polarization axes p = n x L/\n x L\ and q — —n x p. 
The signal received by LISA is then h{t) = y2[F+/i+ + F^'h^]. with: 

F+h+ + Fx/ix = [{paPb - qaqb)h+ + (paQb - qaPb)h^]{iX - t^t\). (3) 
Defining 
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we can write the signal in a standard phase/ amplitude form: 

hit) = J2 cos [27rnft + (;6„p(t) + ^^^(t)] (6) 

where the changing orientation of LISA contributes to a time-dependent phase, 
4>np{t), due to the varying sensitivity to each polarization and the orbital motion 
of LISA about the sun contributes to a time-dependent phase, (pnoit), due to the 
Doppler shift in the received signal. 

The angles 'd and (/? are determined by the source position and orientation by: 
cos-i? = —n ■ L and iamp = {n ■h)/{n ■ a). The polarization phase and Doppler 
phase are: (pnp = tan^^ (— A„/i?„) and (j)nD = 27mf{R/c) sin 6*5 cos (0(t) — (ps) with 
R — lAU, and {Os,4>s) the angular position of the source. The center-of-mass 
trajectory of LISA is 0(t) = 27it/T with T = lyr. 

A second signal is generated from a linear combination of the responses from all 
three arms of LISA which is equivalent to the signal from a detector rotated by 45° 
in the plane of the initial detector [4] . These two signals are used to calculate the 
signal to noise ratio for eccentric binaries. 

The Fourier transform of the signal is h{f) = J2 h^^Kf) with h^''^^ sharply peaked 
about nf. The signal to noise ratio (p) is found by comparing the strength of the 
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signal at each frequency to the expected instrumental noise at that frequency. The 
spectral noise density Sn{f) can be considered nearly constant in the small region 

about nf. Using Parseval's theorem and the fact that the amplitude is slowly 
varying with respect to the harmonic frequency, we have: 

E j:[Sn{nfr' nAl + Bl)dt. (7) 
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We have explored the parameter space of possible eccentric binaries using a 
Monte Carlo approach. The source position is assigned to be the position of a 
globualr cluster selected at random from the galactic globular cluster system. The 
source orientation is assigned with all orientations equally likely. The mass of 
each component of the binary as well as the eccentricity and orbital frequency 
are all chosen with equal likelihood from the following ranges: 8Mq < M < 
18Mq, < e < 1, and O.SmHz < / < O.OlmHz. Out of a total of 300 binaries, 
71 were detectable by LISA with p > 2. Most of these had significantly higher 
values of p. The number of detectable binaries in each bin of orbital period and 
eccentricity are shown in Fig 1. The detection of an eccentric binary relies an 
multiple harmonics appearing in the data stream. The signal strengths from a 
simulation of 100 globular cluster binaries are also shown in Fig 1. 
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FIGURE 1. On the left, a contour plot showing the relative number of detectable binaries for 
a given range of orbital period (in days) and eccentricity. Darker bins contain more detectable 
binaries. On the right, the signal strengths from a simulation of 100 globular cluster binaries. 
Each track represents the harmonics for a single eccentric binary 
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